Enhanced dielectric response by disordered nanoscale/mesoscopic insulators 
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Enhancement of the dielectric response of insulators by disorder is theoretically proposed, where 
the quantum interference of electronic waves through the nanoscale/mesoscopic system and its 
change due to external perturbations control the polarization. In the disordered case with all the 
states being localized, the resonant tunneling, which is topologically protected, plays a crucial role, 
and enhances the dielectric response by a factor 30^40 compared with the pure case. Realization 
of this idea with accessible materials/structures is also discussed. 
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Dielectric property of insulators has been one of the 
most fundamental issues in condensed matter physics. 
Recently, it has been utilized as the ferroelectric ran- 
dom access memory (FeRAM) 0, where thin films of fer- 
roelectrics are integrated onto semiconductor Si wafers. 
The current FeRAM is comprised of a metal-insulator- 
metal structure illustrated in Fig. 1 (a). Electric polar- 
ization is reversed by applying an anti-parallel voltage 
pulse, trasferring from one lead to the other the electric 
charge determined by 2P r with P r being the remnant po- 
larization. By contrast, applying a parallel voltage pulse 
does not yield a polarization change or a charge trans- 
fer. This polarization current is the operation principle 
of the current FeRAM. Especially thin ferroelectric films 
of the nanoscale thickness are used, where the quantum 
effect becomes more important. Therefore, a quantum 
theory of nanoscale/mesoscopic capacitance is called for. 
For applications, it is also often required to achieve (i) 
a magnitude of the polarization larger than 10 /iC/cm 2 , 
(ii) a smaller leak current than 0.1-luA/cm 2 , and (iii) a 
dielectric constant larger than 300 

Generally, the dielectric response is inversely propor- 
tional to the energy gap Eg |2,|3j. Here, Eg is the charge 
excitation gap of the system. Therefore, a small charge 
gap is required to enhance the polarization and the di- 
electric response. One of the most advantageous ways of 
achieving this is to introduce a disorder potential |j] as 
realized in solid-solution systems, since neither dc trans- 
port nor dissipation by leak current occurs due to the 
localization effect 0, 0- I n f ac ^ the localization of the 
wavefunction disconnects a coherent charge displacement 
and undesirably leads to a charge accumulation. This 
difficulty is overcome by taking nanoscale/mesoscopic ca- 
pacitors with the thickness of the order of the localization 
length. 

In this Letter, we develop a quantum theory of 
small-size capacitance, where the coherence of the elec- 
tronic waves through the sample matters. Based 



FIG. 1: (Color online) Ionic dimer system sandwitched by the 
leads, (a) Transmission and reflection coefficients, t (t') and 
r (r) at the left (right) edge of the sample, (b) Ionic level 
alternation (Qi) and bond dimerization (Q2), and effects of 
electric field E on the changes 8Q on Q. 



on this, we propose a disorder-induced quantum- 
mechanical enhancement of the dielectric response in 
nanoscale/mesoscopic ferroelectrics, which provides a 
novel efficient principle of the FeRAM rapidly working 
with low dissipation. It also reveals topological proper- 
ties of the ferroelectrics: in the presence of the disorder, 
a perfect transmission or equivalently a vortex of the re- 
flection coefficients exists as the resonant tunneling in 
the space of the ionic gap and the dimerization. Around 
this vortex center, the most efficient charge transfer can 
be achieved by controling the phase of the reflection co- 
efficient, as has been discussed in the context of quantum 
dots @. 

We consider a one-dimensional (single-channel) insu- 
lating electronic system sandwiched by leads (electrodes) 
as shown in Fig. 1 (a). The extension to higher dimen- 
sional (multi-channel) cases will be discussed later. We 
adopt the Landauer-Biittiker formalism where the 

sample is regarded as a scatterer characterized by the 
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scattering S'-matrix S = 



t 



t* 



The transmittance 



through the sample and the reflectance at the both ends 
are T = \t\ 2 and R = \r\ 2 = \r'\ 2 = 1 - T respectively. 
Then, we employ the Brouwer's formula which has 
been successfully applied to the charge pumping in metal- 
lic quantum dot systems the charge Aq pumped 
from the left lead to the right during an adiabatic change 
of parameters Q along a path C is given by 01 



Aq 



dQ 
2tt 



Im(r*V 



(1) 



Therefore, even in the perfectly reflecting case t = 0, 
charge pumping or transfer occurs from one lead to the 
other by controlling the phase tp of r = \r\e lv [jUT^. 

We take the Hamiltonian schematically shown in 
Fig. [2(b): H = H s + H load , where 

H s = - 2^ o ( c i+l c i + h - c -) 



(2) 



1=1 



for the sample and i?i ca d for the leads and interfaces is ob- 
tained by putting Qi = Q2 — Vi = in H s . Here, ci and 
cj are the annihilation and creation operators of the elec- 
tron at the site i. i n . n . is the transfer integral. Q\ and Q2 
represent the alternations of the local ionic level and the 
bond dimerization, respectively. Vi denotes a local ran- 
dom potential at the site i = 1, ■ ■ • , L with L being the 
number of sites in the sample. The spin degree of freedom 
is omitted for simplicity. H s describes the essence of the 
electric polarization, i.e., the relative displacement of the 
two species of atoms, and describes the ferroelectricity 
in BaTiC>3 0] and quasi-one dimensional ferroelectrics 
such as organic charge transfer compounds TTF-CA 0] 
and (TMTTF) 2 X Experimentally, one can change 

the parameters Q by applying the electric field E along 
the polarization direction and the pressure p. The electric 
field E linearly changes the angle 9 = arctan(Q 2 /<3i): it 
can be viewed that electrons at high and low deinsity sites 
shift in relatively opposite directions, as shown in gray 
arrows of Fig. ^ (b), changing the dimerization Q2 by 
5Q2 oc eEQi. Simultaneously, within each dimer, a level 
difference Q\ changes by 5Q\ oc ~eEQ2 as illustrated 
by black arrows in Fig. ^ (b) . Applying the pressure re- 
duces the magnitude of the gap \J Q\ + Q\ through the 
increase of hybridization. 

Let us start with the pure case = 0). Two upper- 
most panels of Fig. [3 show the reflectivity R = \r\ 2 for 
L = 10001 and the phase tp of r for L — > 00 without ran- 
domness. Here the phase winds by 27r around Q = (0, 0), 
forming a "vortex" at which r — and \t\ = 1. (To 
be precise, the transmittance T behaves as e~ L ^° with 
Co = in. n. /Ego- Therefore the size of the region with 



a large T is of the order of t a . n ,/L for a large L.) To 
understand this result, it is helpful to consider the bulk 
system. Without the random potential the Hamil- 
tonian H s under the periodic boundary condition yields 
a gap Ego = ^\/Q\ + Min{t n n , Q%}, which closes at 
Q = 03]. Namely the vortex corresponds to the gap- 
less case, where the extended state at the Fermi energy 
carries the charge and causes the perfect transmittance. 
This perfect-transmission point Q = governs topolog- 
ical properties of the system in the whole Q plane in a 
non-local way. Along a cycle around Q = 0, the charge e 
is pumped and the polarization changes by ±2ea (a: lat- 
tice constant) according to the Brouwer's formula Eq. 
since t vanishes. This quantized charge pumping |2| is 
completely consistent with the results obtained in the 
periodic system 0, ll^ | by using the Berry-curvature for- 
mulation of the electric polarization 0] . 

Let us introduce the disorder. In the Q plane, the 
region of |Q| < v becomes gapless. The Anderson local- 
ization seriously affects the transport properties 0, 
In one dimension, the effect is pronounced and all the 
states are localized ^3H3;0|- O n the other hand, the 
total vorticity is an integer topological number robust 
under a continuous change of parameters: the phase of 
r (r') should wind by (— )2tt for a large cycle far from 
Q = 0, well within the gapped region There should 
be at least one vortex (r = 0) in the Q plane. Therefore, 
perfect transmittance \t\ — 1 occurs even though all the 
states are strongly localized. 

Figures [3 (a2) and (b2) show our numerical results 
of R and tp for L = 401 with a uniform random dis- 
tribution vi G [— v/2,v/2) for the disorder strength v = 
0.5£ n .n.. The location of the perfect transmittance \t\ = 1 
shifts in the Q space. Besides, the shape of the region 
of relatively high transmittance changes from isotropic 
to highly anisotropic. In this anisotropic "wing" , tp 
changes rapidly. The width Q\y of the anisotropic wing 
decays exponentially as exp(— £/0 with incresing L. 
Here, the localization length £ is obtained from a nu- 
merical calculation of the inverse participation ratio as 
£ = (J2f A) 2 / J2i A 2 : where A is the local density of 
states at the site i. We obtain £ = 76.4 in the case of 
Figs.0(a2) and (b2). In Figs.0(a3) and (b3), we further 
increase the strength of the random potential by taking 
Vi = Si{v + Svi)/2 with v/t n . n . — 2.0 for a smaller sam- 
ple size L — 25, where Sj = ±1 denotes a random sign 
and Sv.i/v £ [—0.05, 0.05] a uniform random distribution. 
This type of randomness mimics the effects of substitu- 
tion in alloys. In this case, the region in the Q plane 
where the gap collapses and vortices distribute expands. 
The localization length £ ^ 3 or 4 is much shorter, and 
the transmittance T is typically ~ 10 -9 , which is practi- 
cally negligible, except at the vortex core. 

Properties around the perfect transmittance can be 
explained in terms of the resonance tunneling 0, 0- 



3 




(a3) 

1.0 



C 



1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 
QiAn.n. XIO " 



v/t M = 2.0 
L-25 



(Ct) 
2.0 

1.5 

1.0 

^ 0.0 

(5-0-5 
-1.0 
-1.5 
-2.0 



1.0 
Q,/tn.„. 



v/t„.„ = 1.6 
L=25 



(b3) 
1.0 



O 



-1.0 . 
2.0 0.0 



(c2) 

2.0 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 
Q,/t„.n. 



0.0 0.2 
XIO" 2 



0.5 1.0 

Q.An.n. 



LI 



-2.0-1.5-1.0-0.5 0.0 0.5 1.0 1.5 2.0 
Q,/t„.n. 



1.5 
1.0 

j| 0.0 

C?-o.s 

-1.0 
-1.5 
2.0 



v/t n .„ .= 2.4 
L = 25 



VI 

i 



n 4J 

1 41 

■ \ » 

: k 

y; 

I — 1 ii 



(d) 



-2.0-1.5-1.0-0.5 0.0 0.5 1.0 1.5 2.0 

Q,/t n .„. 

271 





— ^\ co>(2n+1)7i 


/ S1<S2 


2- — ~" j co<(2n+1)7t 

S1>S2 



FIG. 2: (Color) Reflectance R = \r\ 2 , and the phase ip = 
arg r. (al), (bl): clean bulk system. (a2), (b2): disordered 
system of finite size L — 401 with the random potential of the 
strength v/t n . n = 0.5. (a3), (b3): disordered system with L = 
25 and a stronger randomness v/t n . n . = 2.0 (see the text). 
The color code given in (al)/(bl) also applies to (a2,3)/(b2,3). 
In the white region in (al,2,3), R is the unity within the 
accuracy of 10 -4 . (cl,2): relative dielectric response in the 
presence of the disorder compared with the pure case. The 
disorder strength is v/t n .n. = 1.6 for (cl) and 2.4 for (c2). 
Averages are taken over 10 2 random disorder configurations. 
Inside white bands, there occurs a gradual sign change in the 
dielectric response, (d): Schematic contour mapping of the 
phase cp of r. 



When the Fermi energy coincides with one of the eigenen- 
ergies of localized states, the conductance is resonantly 
enhanced due to tunneling from the usual exponentially 
suppressed value for £ <C L. For an effective potential 
V(x) having two peaks with their heights characterized 
by the dimensionless parameters Si (i = 1,2) 0, 
the reflection coefficient at the left end reads r = 
(tanh Si + tanh S 2 e iu ) / (l + tanh Si tanh S 2 e iu ) , 

with a real constant <f), lu = 2 J , k(x)dx in the semi- 
classical approximation |jj , and an effective wavenumber 
k(x) ~ \jEp — V{x). (Si and S2 are regarded as 
effective parameters when the real potential is much 
more complicated as in our case. They may be directly 
relevant to double-dot systems 22].) Then, a perfect 
transmission \t\ = 1, or equivalently a vortex or r, occurs 
when Si = S2 and to — (2n + 1)tt (n:integer). The latter 
is the Bohr quantization condition that the energy of the 
localized eigenstate is equal to Ep 0- In a cycle going 
around the vortex center in the direction to the arrows, 
the phase of r winds by 2-7T, as shown in Fig. EJd), 
while r' yields the opposite phase winding. Since the 
transmittance is negligibly small except at the vortex 
center, according to the Brouwer's formula Eq. (JTJ, this 
adiabatic cycle pumps nearly one unit charge from the 
left lead to the right, like a bicycle pump 0|. This 
nearly quantized charge pumping was discussed in the 
context of quantum dots [8j . 

When £ <C L, Si and S2 are of the order of L/£ ^> 1, 
and the phases of r and r' abruptly change by 2-7T within 
the range of \u> — (2n + 1)tt| ~ e^ L ^. Along the line 
uj = (2n + l)7r, we obtain r = tanh |Si — S2I, which ap- 
preciably deviates from zero only when Si — S2I > 1, i.e., 
in a small region of the order of A/i n n . ~ £/L. These 
explain our numerical results. In such case, because the 
tunneling rate is exponentially small as Qw ~ in.n.e _i '^, 
an exponentially long time t ~ (fi/tn.n.)e L is required 
to adiabatically pump charge, refleecting the uncertainty 
between energy and time: in order to define the param- 
eters Q with the higher accuracy than Qw-, w e need to 
change Q in a longer time than h/Qw- If the parameters 
are changed in a high frequency u> > Qw /S, the pumped 
charge obtains exponentially decaying factor of Eq /lu as 
a nonadiabatic correctio n I23I. This behavior resembles 
the relaxor ferroelectrics |26ll27l Eif . 

To attain the enhanced dielectric response, the system 
should be located close to the "wing" in the Q space. We 
consider a thin film consisting of many channels, each of 
which is described by the present disordered model but 
with a different profile of random potentials. We can 
design the pattern of the phase <p in the Q plane by the 
disorder. For the alloy model with the random on-site po- 
tentials, vortices are mostly located around Qi = ±v/2 
with the "wing" almost parallel to the Q2 axis. There- 
fore, we can enhance the dielectric response by tuning 
the disorder strength so that the realized Q point is 
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located in the "wing" region. In particular, when T 
is negligibly small, the enhancement factor is given by 
{dip/d9) d [ sordcT / (d(p/d9) pul . c from Eq. QJ, since the elec- 
tric field gives a linaer change in 9 — arctan(<5 1 /<5 2 )- 
Figures HJcl) and (c2) represent the Q dependence of 
the enhancement factor in the case of v/t n , m = 1.6 and 
2.4 for L = 25, respectively, with the number of channels 
being 10 2 . The results reveal that around Qi ~ ±u/2, 
the dielectric response is largely enhanced by a factor 
30 ~ 40 compared with the pure case. Even for the thin 
film with a square shape of a linear dimension larger than 
50 A, which corresponds to N = 10 2 , the disorder-driven 
enhancement of the charge tranfer rate remains robust. 
This allows for a reduction of the applied electric field 
necessary for manipulating the polarization. If we re- 
quire the response time r ~ e L ^/t n n . of the order of 
10 -9 s, we obtain e L ^ < 10 6 with the assumption of 
t n .n. ~ 10 15 s _1 . This also assures a negligibly small 
transmittance T ~ e~ L /? ~ 10~ 6 , and thus the small 
leak current and low dissipation. 

The above scenario can be experimentally verified 
by preparing thin films of solid solution such as 
Pb(Feo.5Nbo.5)03 and Pb(Sco. 5^0.5)03, which contains 
a difference in the local electronic level. Note that the 
solid-solution system having B-site ions in equal ratio be- 
comes ordinary ferroelectrics but not a classical relaxor 
if they are prepared with an adequate slow-anneal pro- 
cess |24j,|25j. It can be regarded as the quantum analogue 
of the relaxor: In relaxor ferroelectrics [ifl EtI Eij | and 
pinned charge density waves [2^, the mechanism is due 
to a classical mesoscopic cluster formation analogous to 
spin glass systems. In addition, the choice of the param- 
eters in the present model are realistic for experiments, 
because the critical thickness for ferroelectricity of Pb- 
based perovskite thin films is in the range of 1.2 ~ 4 
nm [3fJ, |3l| and the capacitor area of a practical device 
(2M-bit FeRAM) is down to 0.4224 /1m 2 j^- 

We briefly mention the relevanec of the present theory 
to the charge pumping through a metallic dot |12J weakly 
connected with two leads. Changing potential heights be- 
tween the dot and the neighboring leads, an unquantizcd 
charge is pumped from one lead to the other [llj. Here, 
the main source of the pumping is a finite transmittance. 
It crosses over to the opposite regime with a vanishing 
transmittance where the pumped charge is nearly quan- 
tized around the resonant transmission 

In conclusion, we studied the quantum-mechanical 
disorder-enhanced dielectric response in nanoscalc insu- 
lators. Phase winding of the reflection coefficient r at 
the vortex exists in the space of ionization and dimeriza- 
tion, and it plays a crucial role in adiabatic charge pump- 
ing and changing the polarization of the system through 
the resonance tunneling, whose position varies with the 
chemical potential. This serves as a new quantum theory 
of the FeRAM or capacitance. 
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